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SU(3) lattice gauge theory is studied by means of an improved action where a 2 x 2 Wilson loop 
is supplemented to the standard plaquette term. By contrast to earlier studies using a tree level 
improvement, the prefactor of the 2x2 Wilson term is determined by minimizing the breaking of 
rotational symmetry detected from the static quark-antiquark potential. On coarse lattices, the 
novel action is superior to the Iwasaki action and comparable with DBW2 action. The scaling 
behavior of the novel action is studied by using the static quark potential and the ratio of the 
deconfinement temperature and the string tension. 

PACS numbers: 11.15.Ha, 12. 38. Aw, 12.38.Gc 



A systematic approach to the non-perturbative regime 
of Yang-Mills theory is provided by computer simulations 
of lattice gauge theory. Regularization is introduced by 
replacing continuous spacetime by a hypercubic lattice 
with lattice spacing a. Improved actions for those lattice 
gauge simulations have attracted much interest in the 
recent past, since they allow for simulations on coarse 
lattices without an overwhelming impact of discretiza- 
tion errors. For example if the infra red regime of QCD 
Green functions is addressed (see e.g. the low lying 

momenta are of order l/Na, where N the number of lat- 
tice points in one direction. For a reasonable amount of 
lattice points, large lattice spacings are highly desirable 
for these purposes. Using the standard Wilson action, 
large lattice spacings give rise to sizable violations of ro- 
tational symmetry. These induce a systematic error to 
the lattice data and severely limit their significance. The 
benefit of improved actions is that they make use of large 
lattice spacings without sacrificing the continuum limit. 

In lattice gauge theory, dynamical degrees of freedom 
are unitary matrices C/^(a;) € SU(3) associated with the 
links of the lattice. The partition function is defined by 



Z 



VU^ exp{-5iatt[f/]} 



In order to recover continuum Yang-Mills theory, one de- 
mands in the limit a — > (adopting minimal Landau 
gauge) 

C/^(x) = exp{zi''4(a;)a} , (1) 
5iatt = E [i^M^^M. + Oi a« + . . .] , (2) 

where are the generators of the SU(3) group, A^^{x) are 
the gauge fields and F^^, is the usual non-Abelian field 
strength tensor. When one designs a realization of 5iatt 



TABLE I: Coupling constants of the "1x2" term contributing 
to the improved action. 





Wilson 


Tadpole improved 


Iwasaki 


DBW2 


/ 





l/20?io 


0.0907 


0.1148 



in terms of the matrices C/^ [x) , it is essential that Sutt 
is gauge invariant for any choice of the lattice spacing 
a, thereby enforcing gauge invariance in the continuum 
limit. The most simple choice for such an action is the 
Wilson action: 
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w',:^\x)} , (3) 
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where WJ^^"^{x) is the n x m rectangular Wilson loop, 
located at x, with extension n in /x and extension m in 
v direction. /3 is related to the bare gauge coupling go 
by /3 = 6/(7q. Universality imphes that any other lattice 
action which meets these requirements also serves the 
purpose in the limit a — s- 0. The crucial point, however, 
is that computer simulations are carried out on finite lat- 
tices consisting of N"^ lattice points. The corresponding 
extension of the lattice in one particular direction is given 
by L = Na. If we insist on L > 1.2 fm, the lattice spacing 
must not be smaller than a — 0.075 fm for a 16* grid. The 
central idea is to construct a lattice action which mini- 
mizes the contribution of the irrelevant terms such as Oi 
in . At the expense of a more complicated lattice ac- 
tion, simulations might be performed on coarse lattices 
and might produce results with little corrections even 
for sizable values a. Several proposals have been made 
for such actions: tadpole improved actions 3], improved 
actions obtained by lattice perturbation theory |^, 
improved actions based on renormalization group stud- 
ies 0,0 1 and "perfect^' actions operating close to the fixed 
point of the theory 
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FIG. 1: The measure 5„ of the breaking of rotational sym- 
metry for several improved actions as function of the lattice 
spacing a in physical units. 



Let us focus on improved actions of the type 
'1. 



^latt 



x.,fj.>y 

1 



f-RetT{2-W'^:^'ix)~W'^^\x)} 



(4) 



Different ideas on improvement result in different propos- 
als for /. More than ten years ago, Iwasaki used pertur- 
bation theory for a renormalization group analysis and 
suggested / = 0.0907. The so-called DBW2 action, de- 
rived from "double blocked Wilson" configurations 0, 
chooses f = 0.1148. The "tadpole improved tree level" 
action p| (see also |T3| for details) favors a (3 dependent 
factor / — 1/20 Mg, where 



-Re {trWl:^\x)) 



1/4 



must be self-consistently determined. See Table for 
a summary. A thorough study of the scaling behav- 
ior of several "1 x 2" actions was recently performed 
by Necco ^l|. It was found at least for small a that 
the DBW2 action shows larger lattice artefacts than the 
Iwasaki action. 



TABLE II; String tension in units of the lattice spacing, aa 
using the novel action with / = 0.056. 



/3 


6.0 


6.2 


6.4 


aa'' 


0.343(3) 


0.247(1) 


0.178(10) 
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FIG. 2: The static quark potential for the novel action (/ : 
0.056) for a 16" lattice. 



based upon the lattice action of type 

'1 
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Retr{l-Wl^\x)] 



f-Rctr{l-Wl::\x) ] 



(5) 



Here, we propose a new type of improvement which is 



Since W^^^{x) possesses the same (lattice) symmetry 
structure than the plaquette, the same symmetry break- 
ing terms show up upon an expansion with respect to 
small lattice spacings. One therefore can hope that for a 
suitable choice of / almost all of the breaking terms can 
be removed. Improved actions containing a W^^^ term 
were considered in |0, where the prefactor / was de- 
termined at tree level, i.e., ftree = 1/64. It turned out 
that already the tree level improved action substantially 
progressed the analysis of the hot gluon plasma. Rather 
than using the tree level improvement, we here pursue 
a non-perturbative improvement: / will be determined 
by minimizing the breaking of rotational symmetry by 
the underlying lattice. Since Oi in (|2Jl explicitly contains 
terms which break rotational symmetry Q , the absence 
of rotational symmetry breaking is directly a measure 
for the suppression of irrelevant terms and, therefore, for 
improvement. 

For an analysis of rotational symmetry breaking, let 
us consider the potential of a static quark-antiquark pair 
separated by a distance vector r. Choosing one of the 
main axis of the cubic lattice, i.e., r = rei, i = 1 . . .3, 
yields the "on-axis" potential V"o„(r). Alternatively, one 
might choose r = r/\/2 (1, 1, 0)"^ (and permutations) and 
r — r/-\/3 (1, 1, l)"^- Contributions from these directions 
to the potential results in the "off-axis" potential. The 
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SU(3) Wilson. b^5.69, 8000 configs 
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FIG. 3: The probability distribution P{y) of the Polyakov line close to Tc for the Wilson action (left). P{y) for the novel 
action employing a 15 x 5 lattice. 



TABLE III: Locating Tc for the novel action 



X Nt 








8^ X 2 


5.718(4) 


0.5381(9) 


0.682(11) 


9^ X 3 


6.21(1) 


0.245(2) 


0.673(5) 


12^ X 4 


6.55(2) 


0.1446(3) 


0.657(9) 


15^ X 5 


6.88(2) 


0.0961(3) 


0.64(1) 



numerical method used here to calculate potentials in- 
volves an overlap enhancement of the Wilson loops with 
the quark ground state along the lines described in 
In fact, a slightly improved version of overlap enhance- 
ment was used, the details of which can be found in 1^. 
Let SV{r) denote the statistical error of the potential at 
distance r. A measure of rotational symmetry violation 
was constructed in 0, i.e.. 

In practice, the "on-axis" potential is well represented by 
the fit 

Von{r) = a + b/r + cr . (7) 

The latter expression provides information on the "on- 
axis" potential at positions r which are not multiples of 
the lattice spacing. Figure shows the measure of rota- 
tion symmetry breaking, 6y, for several actions. Data for 
Sv for the DBW2 and the Wilson case are available in ■ 
My results are in good agreement with these findings. For 
the first time, Sy for the Iwasaki case are reported. We 
find that Iwasaki's choice for / does not completely re- 
move the rotational symmetry breaking clfects at lattice 
spacings as large as ~ 0.1 fm^. On coarse lattices, the 
DBW2 action is less plagued by corrections due to irrel- 
evant terms. By contrast, for < 0.025 fm^ the Iwasaki 
action seems to be less biased by lattice artefacts. My 
findings therefore agree with the results of Necco in . 



Concerning the 2x2 action equipped with tree level fac- 
tor / = 1/64, the improvement on rotational symmetry 
breaking effects is hardly better than in the Wilson case. 

The key idea of non-perturbative improvement is to 
free the parameter / — 1/64 of the 2x2 action, and 
to find a choice for / which is optimal with respect to 
the symmetry breaking effects. One possibility for this 
program is to perform a renormalization group study 
employing lattice perturbation theory, see or using 
blocked configurations, see 0. Here, we adopt a practi- 
cal point of view: For a given (3, we calculate the quantity 
6v for several values /. Starting with the tree level value 
/ — 1/64, we systematically increase / in steps of 10^^, 
and we are monitoring Sy. Once the order a^-terms, the 
terms generating the linear slope of Sy in figure ^ are 
eliminated, we stop the procedure. For this aim, 100 
configurations were generated using a 10^ lattice. Config- 
urations were obtained by virtue of the Cabibbo Marinari 
algorithm [l(| supplemented with micro-canonical reflec- 
tions in order to reduce autocorrelations (details can be 
found in [iSjl). The above procedure for determining / 
is repeated for several values /3 implying that / could in 
principle be a function of (3. However, it turned out that 
the choice 

/ = 0.056(1) (8) 

yields satisfactory results for all (3 values under investiga- 
tions. Note that this value is roughly a factor 3.5 larger 
than the tree level choice, i.e., ftree = 1/64 = 0.015625. 
We find at the end that the "2 x 2" action with / = 0.056 
discussed in the present paper yields results of the same 
quality than the DBW2 action, in particular, in the 
regime of coarse lattices. In addition, the novel action 
seems to be comparable with the Iwasaki action at small 
a. 

Let us finally explore the scaling properties of the novel 
action for a 16^ lattice. For each /3 G {6.0, 6.2, 6.4}, 136 
independent configurations have been generated. The 
static quark-antiquark potential is fitted with the ansatz 
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FIG. 4: Scaling analysis for Tc/^/a for the Wilson action and 
the novel action. 

Q. The results for the string tension are listed in 
table |nj Using these values for aa^, we finally ex- 
press the static potential (as wrell as r) in physical units 
(ct — (440 fm)^ is set for the reference scale). The result 
is shown in figure |21 We observe almost perfect scaling: 
the data points obtained from different /3 values fall on 
top of a single curve. In addition, effects from the break- 
ing of the rotational symmetry by the underlying lattice 
are invisible. 

As a further application, let us study the deconfine- 
ment phase transition at finite temperatures. The order 
parameter of the transition is the expectation value of 
the Polyakov line, i.e.. 




where L is the spatial extent of the lattice, and 
V{x) = \{Uo{x,n) 

n=l 



with Nt being the number of lattice points in Euclidean 
time direction. The temperature T is given by T = 
1/Nta{(3). Since we are only interested in the scaling 
properties of the novel action, we have chosen a physical 
situation which is accessible by a reasonable amount of 
computer time: Rather than to aspire the infinite vol- 
ume limit — > oo, we studied the finite volume case 
L = i/Tc- Thereby, Tc is defined as follows: Let P{y) de- 
note the probability distribution for the variable y in lO . 
At the deconfinement phase transition, the probability 
distribution shows a characteristic double peak structure 
reflecting the first order nature of the transition (see fig- 
ure In the finite volume system, is defined as the 
temperature at which the two peaks of P[y) are of equal 
height. For a given Nt, simulations extracted the criti- 
cal coupling f3c for the novel action and for the Wilson 
action for comparison. The findings for the novel action 
are summarized in table Hill For the largest lattice and 
/3 = 6.91, the auto correlation length was estimated to 
T « 12, whereas 1300 configurations were used to esti- 
mate P{y). Our final results for Td^fo are summarized 
in figure^ 

In conclusions, a thorough study of the non- 
perturbatively improved 2x2 action was performed. 
While the Iwasaki action shows a significant breaking 
of rotational symmetry for lattice spacings larger than 
acP' > 0.25, the present action yields results in this 
regime which are even comparable with those of the 
DBW2 action. On the other hand, scaling violations, de- 
tected from Tcj \fo, are of the same order as those from 
the Iwasaki action. 
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